Flavor decomposition of the elastic nucleon electromagnetic form factors 
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The it- and d-quark contributions to the elastic nucleon electromagnetic form factors have been 
determined using experimental data on G™, G^, G|, and G^. Such a flavor separation of the form 
factors became possible up to 3.4 GeV 2 with recent data on G ™ from Hall A at JLab. At a negative 
four-momentum transfer squared Q 2 above 1 GeV 2 , for both the u- and d-quark components, the 
ratio of the Pauli form factor to the Dirac form factor, F2/F1, was found to be almost constant, 
and for each of F2 and Fi individually, the d-quark portions of both form factors drop continuously 
with increasing Q 2 . 
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Electron-nucleon scattering has been extensively stud- 
ied in two cases. The first case is in elastic scattering 
which is characterized by the electromagnetic form fac- 
tors pp. The second case is in deep inelastic scattering 
(DIS) characterized by the structure functions which ex- 
hibit Bjorken scaling [3J. The study of the proton form 
factors in elastic scattering by Hofstadter et al. provided 
some of the first information on the size of the proton and 
the distribution of charge and magnetization [3J. Deep 
inelastic scattering resulted in the discovery of quarks 
at the Stanford Linear Accelerator Center (SLAC) [I], 
and also taught us about the nucleon's spin structure [5]. 
Elastic and deep inelastic scattering provide complemen- 
tary information: elastic scattering reveals features of 
transverse structure, and DIS provides information re- 
garding the longitudinal momentum distributions of the 
quarks. Taken together, elastic and the deep inelastic 
scattering can both be understood within the broader 
framework of generalized parton distributions which fa- 
cilitate a tomographic picture of the nucleon [BJ. 

Experimental data on the proton Dirac form factor 
F\ [Z] have been found to be in fair agreement with a 
scaling prediction based on perturbative QCD (pQCD), 
Ff cx Q~ 4 |8J, where Q 2 is the negative four-momentum 
transfer squared. However, it has been argued that 
pQCD is not applicable for exclusive processes at ex- 
perimentally accessible values of momentum transfer [9] . 
Indeed, experimental results from Thomas Jefferson Lab- 
oratory (JLab) [TU] for the ratio of the proton Pauli 
form factor F% and the Dirac form factor F[ have been 
found to be in disagreement with the suggested scaling 
F^/Ff cx l/Q 2 .8 J. These same data, however, are in 
reasonable agreement with an updated pQCD prediction 
Q 2 F 2 /F 1 <x ln 2 [Q 2 /A 2 ] [EL] even at modest Q 2 of several 
GeV 2 . Here A is a soft scale parameter related to the size 
of the nucleon. The prediction has the important feature 
that it includes components of the quark wave function 
with nonzero orbital angular momentum. 

In view of these facts it is of significant interest to look 



for the origin of the observed Q 2 -dependence of F$ /Ff . 
We report here on the flavor-separated form factors for 
the up and down quarks up to Q 2 = 3.4 GeV 2 . When 
considering the ratio F 2 /Fi for the d and it-quark con- 
tributions to the nucleon form factors, we find their Q 2 
dependencies to be surprisingly constant. However, when 
combined in the proton form factors, they give the ap- 
pearance of the onset of scaling. It is interesting to note 
that the authors of [llj did not expect the asymptotic 
predictions for the form factors to work at a few GeV 2 
and considered that "... the observed consistency might 
be a sign of precocious scaling as a consequence of deli- 
cate cancellations in the ratio" . 

In the one-photon exchange approximation, the ampli- 
tude for electron-nucleon elastic scattering can be written 
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is the hadronic matrix element of the electromagnetic 
current operators for the proton (neutron) . Here, the ne- 
glect of heavier quarks in this context is supported by ex- 
perimental data on parity non-conserving polarized elec- 
tron scattering from the proton [T2J. While we can not 
evaluate the matrix elements of uj^u and d/y^d explicitly, 
from symmetry considerations we know that the matrix 
element shown in Eq. [I] must have the form (considering 
the proton for definiteness) 
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where p(k) and p(k') are the proton Dirac spinors for the 
initial and final momenta k and k', respectively. The def- 
inition of the neutron form factors F"(Q 2 ) and -F 2 "(Q 2 ) 
follows similarly. 

The JLab data for G p /G p from Refs. [101 were used to 
plot Sp = Q 2 F P /Ff in the upper panel of Fig. Ill which 
also shows the prediction [TT] at A = 300 MeV. Data 



2 




4 5 6 

Q 2 [GeV 2 ] 

FIG. 1: The ratio of the Pauli and Dirac form factors, multi- 
plied by Q 2 , S =Q 2 Fi/F2, vs. the negative four-momentum 
transfer squared Q 2 . The upper panel shows S p for the proton 
and S„ for the neutron using data from Refs. [13-18], as well as 
the curves of the prediction [11]: ln 2 [Q 2 /A 2 ] for A=300 MeV 

The bottom 
and Sd for the 



which is normalized to the data at 2.5 GeV 2 
panel shows the individual flavor quantities S u 
u and d quarks, respectively. 



on G™/G™ for the neutron up to Q 2 =3.4 GeV 2 were re- 
cently published by Riordan et al. [13] . For the first time, 
it is possible to examine the behavior of the neutron ratio 
Fz/F™ in the same Q 2 range as that where the interest- 
ing behavior was first seen for the proton |10) . Using the 
data of Riordan et al. as well as those of Refs. [14-18], we 
also show in Fig. [ljthe quantity S n =Q 2 F?/Ff. Scaling 
of S n is clearly not evident at the lower Q 2 values shown, 
although the data do not rule out this type of behavior 
at a moderately higher Q 2 



Thus far, by discussing F± n ' and F^" 1 ' we are ex- 
plicitly examining the behavior of the matrix element of 
the electromagnetic operators (|u7 m m + ^dry^d) in the 
proton (neutron). If we assume charge symmetry (thus 
implying (p|u7^ii|p) = (n|d7 M d|n)), it is possible to per- 
form a flavor decomposition of the form factors F^™- 1 
and F2^ n \ and construct form factors corresponding to 
the matrix elements of uj^u and d'jn d individually |19j . 
Here we use the relations 
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In what follows, we use the convention that F^ 2 ) and 
-F\ d ( 2 ) TC ^ CI to the up and down quark contributions to 
the Dirac (Pauli) form factors of the proton. At Q 2 =0, 



the normalizations of the Dirac form factors are given by: 
F"(0) = 2 (Ff(0) = 1) so as to yield the normalization 
of 2 (1) for the u (d)-quark distributions in the proton. 
The normalizations of the Pauli form factors at Q 2 =0 are 



given by F%(0) 



x q j 



where k u and Kd can be expressed 



in terms of the proton (n p ) and neutron (n n ) anomalous 
magnetic moments as 



2^r) 



+1.67 and Kd 



2k, 



-2.03. 



Having defined the flavor-separated Dirac and Pauli 
form factors, we can also define the quantities 

S u = Q 2 F?/F? and S d = Q 2 F 2 d /F?, 

which we have plotted in the bottom panel of Fig. [I] Each 
individual data point corresponds to an experimental re- 
sult on G^/G^j from Refs. [13-18]. Only the uncertainties 
in the ratio G ™ IG n are included in the error bars of the 
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flavor-separated results because the other form factors 
(calculated with the Kelly fit (20]) are known to much 
higher accuracy, albeit dependent on the particular pa- 
rameterization chosen. The behavior we see is completely 
different from that of the proton and the neutron. There 
is a striking lack of saturation, and indeed the variation 
of S u and Sd with Q 2 appears to be quite linear. It is in- 
teresting also that the slope associated with the d quark 
is about six times larger than that of the u quark. When 
we consider the matrix elements of u^^u and d'y^d indi- 
vidually, the relationship between the Pauli and the Dirac 
amplitudes is quite different from when we consider the 
sum of the amplitudes that results in the full hadronic 
matrix element (Eq. 2). 

While it is instructive to plot S u and Sd so that we can 
compare them directly with the widely discussed S p for 
the proton, the inclusion of the factor of Q 2 masks the 
detailed behavior as Q 2 approaches zero. We thus plot 
in the top two panels of Fig. [2] the quantities k" 1 F% /F™ 
and k^F^/F^. Here, a second aspect of the behav- 
ior of the flavor decomposed form factors appears that is 
quite intriguing. These ratios are relatively constant for 
Q 2 greater than ~ 1 GeV 2 , but have a more complex be- 
havior for lower values of Q 2 . This might be interpreted 
as a transition between a region where the virtual pho- 
ton coupling to the three-quark component in the wave 
function dominates (higher Q 2 ) and a region where the 
inclusion of a coupling to a five-quark component is es- 
sential (lower Q 2 ). We note also that the ratio F2/F1 
for the proton does not show a different behavior above 
and below 1 GeV 2 as one can see in the bottom panel of 
Fig. [2j The calculation of the form factors in a relativis- 
tic constituent quark model (RCQM) [3T] (shown by the 
blue curves in Fig. [2]) deviates considerably from the data 
which illustrates the discriminating power of the flavor 
separated form factors. The empirical Kelly fit (which 
predates Ref. [T3]), corresponds to the black curves, and 
is in reasonable agreement with the data, particularly at 
lower Q 2 . 
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FIG. 2: The ratios k' 1 f£ / 'F? , K^Fg/F? and K^Fi/Ff vs. 
momentum transfer Q 2 . The data and curves are described 
in the text. 

The form factors Ff, Ff, and F 2 " are shown in 
Fig. [3j all multiplied by Q A for better clarity in the high- 
er 2 range. The values are given in Table F 

TABLE I: The flavor contributions to the proton form factors, 
obtained using G"/G™ form factor data from Refs. [13-18] and 
the Kelly fit [2D] for the other form factors. The Q 2 values 
are given in GeV 2 . 



Q 2 


Ref. 


Fx 


pd 


F% 


Fi 


0.30 


[17] 


1.075(6) 


0.505(12) 


0.716(6) 


-0.995(12) 


0.45 


PS 


0.853(6) 


0.377(12) 


0.515(6) 


-0.777(12) 


0.50 


p] 


0.789(6) 


0.332(12) 


0.473(6) 


-0.708(12) 


0.50 


PI 


0.789(4) 


0.340(7) 


0.463(4) 


-0.713(7) 


0.59 


na 


0.695(6) 


0.283(13) 


0.394(6) 


-0.617(13) 


0.67 


m 


0.628(6) 


0.249(12) 


0.342(6) 


-0.552(12) 


0.79 


m 


0.544(8) 


0.206(15) 


0.283(8) 


-0.467(15) 


1.00 


pi 


0.434(5) 


0.154(10) 


0.211(5) 


-0.357(10) 


1.13 


|18| 


0.379(3) 


0.124(5) 


0.183(3) 


-0.298(5) 


1.45 


[18] 


0.290(3) 


0.093(6) 


0.128(3) 


-0.213(6) 


1.72 


|13| 


0.2257(22) 0.0529(43) 0.1103(22) 


-0.1429(43) 


2.48 


|13| 


0.1380(18) 0.0278(35) 0.0632(18) 


-0.0707(35) 


3.41 


|13| 


0.0851(12) 0.0131(24) 0.0370(12) 


-0.0337(24) 



Up to Q 2 « 1 GeV 2 there is a constant scaling fac- 
tor of ~2.5 for Fi and ~0.75 for F2, between the u- and 
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FIG. 3: The Q 2 -dependence for the u- and d-contributions to 
the proton form factors (multiplied by Q 4 ). The data points 
are explained in the text. 

(i-quark contributions. Above 1 GeV 2 the <i-quark con- 
tributions to both nuclcon form factors multiplied by Q A 
become constant in contrast to the u-quark contributions 
which continue to rise. These experimental results are in 
qualitative agreement with the predictions for the mo- 
ments of the generalized parton distributions reported in 
Ref. [22] • It is interesting to note that the ^-contributions 
correspond to the flavor that is represented singly in the 
proton, whereas the w-contributions correspond to the 
flavor for which there are two quarks. In the framework 
of Dyson-Schwinger equation calculations, the reduction 
of the ratios Ff/F" and F 2 d /F 2 n at high Q 2 is related to 
diquark degrees of freedom [23] • The reduction of these 
ratios has the immediate consequence that S p has its ob- 
served shape despite the fact that S u and Sd are almost 
linear with Q 2 . 

Another representation of the Dirac form factor is the 
infinite momentum frame density, p D , given by the ex- 
pression p D {b) = J(QdQ/2ir)J \Qb)F 1 {Q 2 ) [21], where 
J is the zeroth order Bessel function and b is the im- 
pact parameter. The faster drop off of the d-quark form 
factors in Fig. [3] implies that the u quarks have a signif- 
icantly tighter distribution than the d quarks in impact- 
parameter space, as was noticed in Ref. [25] . 

In summary, we have performed a flavor separation 
of the elastic electromagnetic form factors of the nu- 
cleoli. We find that for large Q 2 the d-quark contri- 
butions to both proton form factors are reduced rela- 
tive to the it-quark contributions. We find also that the 
Q 2 -dependencies of the flavor-decomposed quantities S u 
and Sd are relatively linear in contrast to the more com- 
plicated behavior of S p and S n . This linearity is due 
to the fact, as yet unexplained, that the ratios F^/F^ 
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and F^/Ff are constant within experimental errors for 
Q 2 > 1 GeV 2 . At Q 2 < 1 GeV 2 , however, these same 
ratios show significant variation. Given the linearity of 
S u and Sd, it is quite clear that the precocious scaling of 
the proton form factors and the consistency of the pro- 
ton data with the updated pQCD description of Ref. [TT] 
are the result of the different behaviors of the u- and d- 
quark contributions to the proton form factors. Further 
measurements of G™/G™ j^S] will allow us to extend the 
flavor decomposition to Q 2 =10 GeV 2 and to explore the 
Q 2 range over which the apparent constant behavior of 
F$/F? and Fg/F? persists. 
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